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__ ' We study the pion- nucleon system in s-wave in the framework of lattice 

QCD in order to gain new information on the nucleon excited states. We 
perform simulations for Uf = 2 mass degenerate light quarks at a pion mass 
of 266 MeV. The results show that including the two-particle states drasti- 

j^ i cally changes the energy levels. The variational analysis and the distillation 

approach play an important role in the extraction of the energy levels. The 

Q_|, phase shift analysis allows to extract information on the resonance nature 

D ■ of the observed states. 

(N 
> 



O 



X 



PACS numbers: 11.15.Ha, 12.38.Gc 



1. Introduction 



Almost all the hadrons that constitute the QCD spectrum are unstable 
under strong interactions. Lattice QCD calculation have been tradition- 
ally treating these states as stable states, without taking into account their 
^T) • resonant nature. Only recent studies have made exploratory steps in this 

direction, successfully studying mesonic resonances [Il[2l[3llH[5l[6l[7l[8]. 

We study for the first time the coupled pion-nucleon system explicitly 
including the two particles in our simulations [^ . This work is motivated by 
the fact that lattice hadron spectroscopy does not satisfactory reproduce the 
a ■ negative parity sector of the nucleon states. The physical spectrum consists 

of two resonances A^*(1535) and A^*(1650). So far lattice simulations ^Q\ 
nn [T2l [T3l [T^ have measured in this channel two low-lying states that are 
assigned to the two resonances, even though the lower measured state lies 
below the physical value of A^*(1535) |15] . All these simulations considered 
only 3-quark interpolators that should in principle couple to meson-baryon 
states via dynamical quark loops. However this coupling seems to be weak 
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and meson-baryon interpolators have to be explicitly included in the set of 
operators in order to achieve a complete study of these resonances. 

The negative parity resonances of the nucleon couple to A^vr in s-wave, 
but it is not the only decay channel: Af*(1535) -^ Nr], iV*(1650) -^ iV??, KK 
[16j . However, we work at an unphysical pion mass (m-Tr = 266 MeV) that 
prevents this channel from being in the influence region of the two reso- 
nances. We therefore simulate on our lattice the coupled channel of a 3- 
quark nucleon together with the (4+l)-quark iVvr system in the rest frame. 
The results presented here have already been published in [S]. 

2. Tools and setup 

2.1. Variational analysis 

The energy levels of the nucleon and the A^vr system are determined 
using the variational method [T71 [HI [19] . We measure the Euclidean cross- 
correlation matrix C{t) between different interpolators Oi{t) 

aj{t) = (O,(t)O,(0)) = Y,iOi{t)\n)e-^-\n\O,{0)) , (1) 

n 

and then solve the generalized eigenvalue problem 

C{t)Un{t) = Xn{t)C{to)Un{t) (2) 

to disentangle the eigenstates with the eigenvalues A„(i,to) ~ Q-^nit-to)_ 
The energy values of the eigenstates are determined by exponential fits. 
The fit range is indicated by a plateau-like behavior of the effective masses 

£;„(t) = log[A„(t)/A„(t + l)]. 

2.2. Distillation 

The evaluation of the correlation matrix in the case of (4+1) quarks 
turns out to be hardly accessible to traditional techniques, due to the large 
amount of different diagrams involved. The distillation method |20] allows 
to evaluate partially disconnected diagrams within an affordable amount of 
computer time. The quark sources are smeared using a truncated expansion 
of the 3D Laplacian operator 

q{x) — )■ S{x, X )q{x ) = N^ v^{x)v^ (x ) q{x ) . (3) 

The correlation function for the 3-quark nucleon operator reads 

C(i,„fe,Wc) = <^w(i,i,A:) T{i,i')T{j,j')T{k,k') 4^Ji',j',k') (4) 
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where the perambulators r(?i, m) are quark propagators from source eigen- 
vector f™ to sink u". The functions (p include all the information on the 
Dirac structure of the specific interpolator 

<Pt,„kihJ^k) = ^eabcDvl{x)vl{x)v^ix), (5) 

X 

where D carries all the Dirac indices. 



2.3. Interpolators 

The set of interpolators has to be as complete as possible in order to 
reliably extract the spectrum. We use the nucleon interpolator 

iV» (p = o) = Y^ eabc P± r? Ua (x ) ul (x) F^'^ dc (x) , (6) 

X 

where (ri,r2) = (1,(775), (75,0), (il, (77475) and each quark source is 
smeared combining N^ = 32 and 64 eigenvectors. For the A^vr system we 
use 

Ar7r(p = 0)=75iV+(p = 0)7r(p = 0), (7) 

and we project to isospin 1/2: Onw = P"^^ + V^nir^ with 

7r°(0) = -^ '^{Uaix)-f5Ua{x)-daix)^5da{x)} , 7r+(0) = ^ fia(^)75^ia(^) • 

X X 

(8) 

2.4. Interpretation of the energy levels 

Once the energy levels are computed, one can relate the measured spec- 
trum of the coupled system to the physical resonances. In the elastic region 
Liischer's formula gives a relation between the discrete energy levels (of 
the rest frame system, for a discussion of meson-baryon systems in moving 
frames see [21] ) and the phase shift in the continuum [181 [22] , 

3/2 
tan%)= ^""f \\ ^ (9) 

where the generalized zeta function Zi^ is given in [22] and 

*L ^2 [s - {rriN + m^f][s - {niN - m^)^] 

'i = P^^P= Ts ('°) 

with s = (En)'^. Given some phase shift model, eq. (|9|) can be numerically 
inverted to obtain the expected energy levels for the two interacting particle 
system (Fig. [S] rhs). For a different method to predict the expected energy 
levels in finite volume see e.g. 
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3. Results 

We use 280 configurations generated for two flavors of mass-degenerate 
light quarks and a tree level improved Wilson-Clover action, m.,^ = 266 
MeV, a = 0.12 fm, 1/ = 16^ x 32 



3.1. One particle sector 

Using a set of 3-quark interpolators we reproduce the usual observed 
spectrum [15]. In the positive sector we observe the nucleon ground state 
at jnjv = 1068(6) MeV and another state that lies far above the physical 
Roper. In the negative sector we observe two nearby levels, the lowest lying 
below 7V*(1535) (Fig. [H Ihs). 
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Fig. 1. Effective mass values for A^_. Left: Results for 3 quark interpolators. The 
blue line denotes the non-interacting Ntt state. Right: Results for the N^ and 
N+TT coupled system. 



3.2. Coupled N and Ntt system 

First we compute the energy level for the two particles propagating 
independently (i.e., the threshold) and we find that it is overlapping with 
the first of the two levels measured in the single-particle approach (Fig. [H 
Ihs). When On-k is included a new energy level appears and the effective 
energy levels of the A^vr system show less fluctuations compared to the 3- 
quark case. The lowest level now lies slightly below the A^vr threshold, 
a feature typical for attractive s-wave and a finite volume artifact. The 
next-higher two levels now lie approximately 130 MeV above the physical 
resonance positions of A^*(1535) and A'^*(1650), similar to the shift of the 
nucleon ground state for this value of m-,^. 
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3.3. Phase shift analysis 

A comparison with the expected energy levels obtained inverting Liischer 
formula Q and assuming a single elastic resonance parameterization, show 
excellent agreement (Fig. [2l rhs). Assuming a Breit-Wigner shape for the 
first resonance we can also extract the resonance mass: mn = 1.678(99) 
GeV. 
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Fig. 2. Lhs: Comparison between the experimental masses of the negative parity 
nucleon resonances, the energy spectrum obtained from the single particle analysis 
and the results from the coupled N and Ntt system. Rhs: Energy levels expected 
for the interacting N-k system obtained inverting the Liischer formula and assuming 
a Breit-Wigner parametrization for 7V(1530). 



4. Conclusions 

This study is intended to shed some light on the excited energy levels 
of the nucleon spectrum, which still represents an outstanding challenge for 
lattice QCD. We find that including meson-baryon interpolators is indeed 
needed for a reliable picture of the N- spectrum. The study of two particle 
systems on the lattice improves our understanding of LQCD and this work 
is a first step into that direction. 
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